Real time electron dynamics in an interacting vibronic molecular quantum dot 
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We employ the time-dependent non-crossing approximation to investigate the joint effect of strong 
electron-electron and electron-phonon interaction on the instantaneous conductance of a single 
molecule transistor which is abruply moved into the Kondo regime by means of a gate voltage. 
We find that the instantaneous conductance exhibits decaying sinusoidal oscillations in the long 
timescale for infinitesimal bias. Ambient temperature and electron-phonon coupling strength influ- 
ence the amplitude of these oscillations. The frequency of oscillations is found to be equal to the 
phonon frequency. We argue that the origin of these oscillations can be attributed to the interference 
between the emerging Kondo resonance and its phonon sidebands. We discuss the effect of finite 
bias on these oscillations. 
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I. INTRODUCTION 

Molecular transport junctions that consist of a 
molecule inserted between contacts long fascinated both 
physicists and chemists since it has been proposed more 
than three decades ago that they could be used one day as 
building blocks for future electronic devices. The emer- 
gence of the field of molecular electronics originates from 
this suggestion^. The nanotechnology revolution which 
took place in mid 1980's paved the way for precise con- 
trol on nanostructures and enabled to perform the exper- 
iments that were beyond anyone's realm previously^. Ex- 
perimental confirmation of the Kondo effe ct, which was 
predicted two decades ago theoretically^^, first on quan- 
tum dotiP^ and then on single molecule transistors^! 
opened an avenue for the possible merger of the fields 
of spintronics and molecular electronics^. The goal of 
the field of molecular spintronics is to investigate spin- 
dependent transport in molecular electronic devices^. 
Therefore, the Kondo effect acts as a supplemental spin- 
dependent mechanism that might strongly affect the spin 
transport in quantum dots. 

In his groundbreaking worlP^l, Jun Kondo discovered 
that resistivities of bulk metals which contain magnetic 
impurities with localized unpaired spins would be en- 
hanced at low temperatures. This enhancement was 
dubbed as the Kondo effect later on. The underlying mi- 
croscopic mechanism for the Kondo effect has been iden- 
tified as the formation of a spin singlet resulting from the 
interaction of the unpaired localized electron and contin- 
uum electrons near the Fermi level in the metaP^. 

Quantum dots are artificial atoms that can contain an 
integer number of electrons^. It is imperative to confine 
odd number of electrons in a quantum dot to observe the 
Kondo effect. This gives rise to a net spin in the dot just 
like in a bulk metal. However, coupling of this net spin 
to the fermionic bath in the metallic leads opens a new 
transport channel in the quantum dot and conductance 
is enhanced at low temperatures instead of resistivity in 



bulk metals. A sharp many-body resonance at the dot 
density of states pinned to the Fermi level of the contacts 
is responsible for the enhancement. This conductance 
enhancement is a robust way of obtaining current flow 
in odd number Coulomb blockade valleys and thus to 
prevent current inhibition in single electron devices. 

Sudden shifting of the gate or source-drain voltage has 
been studied in detail previousl y 15 " 18 ! and three unique 
timescales have been clearly identified in the ensuing 
transient currentS^HUl. Non-Kondo timescale is associ- 
ated with the initial fast rise of the current accompa- 
nied by reshaping of the Breit-Wigner resonance. On the 
other hand, formation of the Kondo resonance and its 
reaching a broad quasismooth structure with a linewidth 
on the order of Kondo temperature is the hallmark of 
longer Kondo timescale. Splitting of the Kondo reso- 
nance for finite bias corresponds to the third and longest 
timescale. This timescale can be determined by measur- 
ing the decay rate of the split Kondo peak oscillations^. 
Later studies focused on the asymmetric coupling of the 
dot to the contacts and concluded that an interference 
between the Kondo resonance and the van Hove singular- 
ities in the density of states of the leads may giv e rise to 
sinusoidal oscillations in the transient currenlP^. Recent 
extension of the diagrammatic Monte Carlo method^ to 
nonequilibrium impurity problems^ 5 -! confirmed the sen- 
sitive dependency of the transient current on the band- 
width of the contacts^. 

The history of electron-phonon interaction in molec- 
ular transport junctions is quite long as it has been 
summarized in a recent review^. Phonon assisted elec- 
tron transport in molecular transport junctions can be 
broadly classified based on the relative time and en- 
ergy scales in the process. The strength of the electron- 
phonon interaction is judged relative to the molecule- 
electrode coupling. In this respect, weak and strong 
electron-phonon coupling regimes emerge. 

The former regime, namely the weak electron-phonon 
interaction, corresponds to nonresonant phonon assisted 
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electron tunneling enco untered in inelastic electron trans- 
port spectroscop y 28 ^ 29 !. The development of scanning 
tunneling microscope and spectroscopy proved to be 
an invaluable tool for inelastic electron transport spec- 
troscopy to define and characterize the conductance 
properties of molecular specie s. It is justified to use 
Migdal-Eliashberg theory-22ED i n this regime. As a 
result, second order perturbation theory on electron- 
phonon coupling over the Keldysh contour leads to Born 
approximation. Various versions of the self-consistent 
Born approximation has been used in several theoreti- 
cal studies^MH. 

The latter regime corresponds to resonant tunnel- 
ing which involves longer electron lifetime and stronger 
electron-phonon interaction. Perturbation theory fails 
in this regime and polaron is formed in the junction. 
Signatures of the resonant tunneling are manifested as 
sidebands of the main Kondo resonance in differential 
conductance^. There are several studies investigat- 
ing this case in steady state. Preliminary accounts of 
the electron-phonon coupling computed the dot Green's 
function using the equation of motion method^ and 
treating the contacts as unaffected by the bosons in 
the wide band limitP^. They both concluded that the 
electron-phonon interaction results in sidebands only in 
one side of the main elastic peak in density of states. In 
later studies, it has been suggested that this is due to 
an invalid approximation employed in calculating the re- 
tarded Green's function and the phonon sidebands should 
appear in both sides of the elastic pealpSHSU Pertur- 
bative renormalization group calculations confirmed this 
latter conclusion^]. 

In strong electron-phonon coupling coupling regime, 
the approach invoking Lang-Firsov canonical transfor- 
mation and non-crossing approximation also found side- 
bands in both sides of the main Kondo resonance 
in time averaged ac conductance, however this last 
method concluded that the zero bias Kondo reso- 
nance has been greatly suppressed^. Recently, inves- 
tigation of the electron-phonon interaction in Kondo 
regime using nonequilibrium equation of motion method 
has found that increasing electron-phonon interaction 
strength gradually destroys the Kondo effecfW This has 
been attributed to the destruction of the coherence in the 
system by the electron-phonon interaction and the shift 
of the energy level due to rearrangement of the phonons. 

In this paper, we will investigate a scenario in which a 
molecular quantum dot is suddenly shifted from a posi- 
tion well below the Fermi level of the leads to a position 
where the Kondo effect is present by means of a gate 
voltage. Previous studies unambiguously demonstated 
that the strong electron-electron interactions give rise to 
a Kondo resonance pinned to the Fermi level of the con- 
tacts decorated with sidebands on each side due to the 
electron-phonon coupling. However, little is known about 
the real time electron dynamics of this system as a re- 
sponse to abrupt perturbations except a recent study us- 
ing the mean field approximation for a spinless modeP^. 




FIG. 1: This figure depicts the density of states of the left 
and right contacts alongside with the density of states of the 
quantum dot in the final state schematically. 



The set-up of the system under consideration is shown 
schematically in Fig. [I] The goal of this paper is to fill 
this gap by reporting the instantaneous conductance af- 
ter an interacting vibronic dot's energy level has been 
moved to its final position. 



II. THEORY 

We model this device by a single spin degenerate level 
of energy e^ot attached to leads through tunnel barriers 
and coupled to a single phonon mode. This model is re- 
ferred as Holstein Hamiltonian. In this paper, we will 
be concerned with the strong electron-phonon coupling 
regime where the electron-phonon term can be discarded 
with a canonical transformation. We carry out the auxil- 
iary boson transformation for the resulting Hamiltonian 
where the ordinary electron operator on the dot is rewrit- 
ten in terms of a massless boson operator and a pseud- 
ofcrmion operator. U — > oo limit is obtained by imposing 
the condition that the sum of the number of bosons and 
the pseudofermions is equal to unity. 

The aforementioned Hamiltonian has three pieces de- 
scribing the contacts, quantum dot and the tunneling 
process between them and it can be written as 

H(t)=H c + H D {t) + H T (t) (1) 

where 

He = y^(efcq — Ha)<J kaa Ckuo 

H D {t) = Y} e dot{t) + \{a + a})]dld a + Ud\d^d\d i 

a 

-\-oj$a/ a 

H T (t) = J2(Vka(t)cl aa d a + h.c). (2) 
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In this Hamiltonian, dj r (d cr ) and c \ aa {cka<j) with 
a=L,R create(annihilate) an electron of spin a in the 
quantum dot and in the left(L) and right(R) leads re- 
spectively. Vk a and [i a arc the tunneling amplitudes and 
chemical potentials for the left and the right leads. a^(a) 
creates(annihilates) a phonon within the quantum dot. A 
is the strength of the electron-phonon coupling and uj is 
the phonon frequency. In this paper, we will use atomic 
units with h = fee = e = 1. 

We will take the hopping matrix elements to be equal 
with no explicit time and energy dependence. Under this 
assumption, broadening of the dot level can be parame- 
terized as r(e) = Tp(e) where T is a constant given by 
f = 27r|F(e/)| 2 and p(e) is the density of states function 
of the contacts. We will use parabolic density of states 
in both contacts with same bandwidth. 

If the electron-phonon coupling is weak, electron- 
phonon coupling term can be treated perturbatively. 
In this paper, we will be concerned with the oppo- 
site case, where electron-phonon coupling is sufficiently 
strong compared to the tunnel couplings. In this regime, 
perturbative solution fails and a suitable canonical trans- 
formation must be applied to eliminate the electron- 
phonon coupling term. The most versatile choice is the 
unitary Lang-Firsov canonical transformation^ given by 



S = exp 



A 

W 



(3) 



This transformation gives 



A 



SaS^ = a d\d a 

Sd a S^ = d a X 
where the operator X is given by 



X = exp 



— ( t_ ) 

UJQ 



(4) 



(5) 



The electron operators d a (d\) in Hx(t) are multiplied 
with X(X') as a result of this transformation indicating 
that the tunneling electrons create and destroy a phonon 
cloud. This leads to polaron formation at the junction. 

Under this canonical transformation, dot Hamiltonian 
turns into 



H D (t) = SH D (t)S^ 



X 2 



J2 \^ot(t)- ^ ) s a d a 



u ■ 



2A 2 

UJQ 



d^dfd^d^ + uj a a) a. 



(6) 



It is clear from this result that the dot level and the Hub- 
bard interaction strengths are renormalized as Cdot{t) = 
^dot(t) — (A 2 /w ) and U — U — (2A 2 /w ) respectively as 
a result of this transformation. 



In the strong electron-phonon coupling regime, i.e. 
when a polaron is formed at the junction, mean-field the- 
ory can be applied and the expectation value of operator 
X given by 



< X >= exp 



- )2 / 




--a [N ph .- 





(7) 



can serve as a substitute for the operator. 

Even though U has been renormalized to U = U — 
(2X 2 /uj ), it is still positive and it overwhelms the line 
width r for realistic systems. Hence, typically one takes 
U — >• co. This forbids double occupancy of the dot level. 
The downside of this advantage is that the standard dia- 
grammatic techniques are not applicable anymore. This 
problem can be circumvented by introducing a mass- 
less boson operator and pseudofermion operator on the 
molecule. The original electron operators can be written 
in terms of these operators as 



d„{t) - tf(t)f a (t) 

dt(t) = ft{t)b{t) 



subject to the requirement 



^ "j fafa 



(8) 



(9) 



which ensures the single occupancy of the dot level. The 
resulting slave boson Hamiltonian becomes 

H(t) = Y(e ka - fi a )cl atT Cka<j + tdotflfa + uWa 



+£(vui)cL CT /.& t + ^.), 



(10) 



where the tunnel coupling and dot level are renormalized 
as V ka (t) = V ka (t) < X > and Z dot {t) = e dot (t) - (A 2 /w ) 
respectively. 

We then invoke the well established non-crossing ap- 
proximation(NCA) to determine the pseudofermion and 
slave boson self-energies. NCA has been shown to give 
accurate results for dynamical quantities save for temper- 
atures below T/Tk ~0.1 or finite magnetic fields. These 
two problematic cases will not be dwelled on in this pa- 
per. 

The net current flowing through the device can be 
calculated from the resulting Green's functions. We 
will write the net current as I(t) — — In(t), 

where lL(i)(In(t)) represents the net current from the 
left (right) contact through the left (right) barrier to the 
quantum dot. The most general expression for the net 
currenp3 has previously been derived by using pseud- 
ofermion and slave boson Green's functions^. Here, we 
will adapt it to the present situation. For symmetric 
coupling, the net current is reduced to 

I{t) = -2tlm{ I dhG^t^ifdt-^)- f R (t-h))) 

(11) 
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and the retarded Green function can be expressed as 

G R (t,h) = -i8(t-t 1 )[<d ( ,(t)di(t 1 )><X(t)Xi(t 1 )> 
+ < d£(iiK(t) >< XtfaWt) >]. (12) 

It turns out that the correlators of operators X(t) and 
X^(ti) can be evaluated precisely. They are given by < 
X(t)Xi(h) >= e^* 1 "*) and < X^t!)X(t) >= e^*^ . 
In these expressions, the phase factor is given by 

0(ti-i) - -5^.(1 - e — o(ti-O) + 

(iV^. + ^Cl-e^^-*))], (13) 

where g is defined as g — ^ and N v h , given by Bosc- 
Einstein distribution N P h, = wsk — function, represents 

e k B T _ 1 

the average number of phonons for temperature T and 
phonon frequency ujq. It has been pointed out before 
that the correlators are approximately equal only at suf- 
ficiently high temperatures such that N p h. is much larger 
than unity. At low temperatures, taking them equal may 
produce erroneous results such as producing phonon side- 
bands only on one side of the main Kondo peak. There- 
fore, in this paper we will not resort to such an approxi- 
mation and keep the correlators in separate terms. 
The retarded Green's function can then be written as 

G R {t,h) = -i6(t-t 1 )[<d„(t)4(t 1 )>e^-V 

+ < 4 (*i)d«r(*) > e*'-* ] 
= -ie{t-t l )[G> (ial {t,t l )et'^- i) 

+G<Jt,t 1 )eW-^]. (14) 

Using the slave boson and pseudofermion decomposition 
of the original fermion operators on the molecule, the 
retarded Green's function can be recasted as 

G R (t,h) = -id(t-t 1 )[G R eudo (t,t 1 )B<(t u t)e'l'^ 
+G<e udo (t,h)B R (t 1> t)e^-^] (15) 

The real-time coupled integro-differential Dyson equa- 
tions for the retarded and less than Green's functions 
are computed in a cartesian two-dimensional grid. The 
values are stored in a matrix and the matrix is prop- 
agated diagonally in time. The phase factors remain 
attached to the pseudofermion Green's functions during 
this procedure in order to incorporate the phonon effects 
properly^. The instantaneous conductance is obtained 
by performing an integration over the lowest row of the 
matrix using the expression 

I(t) = 2TRe([ dt 1 (G^ e „ do (i,t 1 ) J B<(t 1 ,t)e^- t ) 

+G< eud0 (t,h)B R (t 1 ,t)e^-^) x 

(/i(t-*i)-M*-ti))), (16) 

where /^(t — t\) and f^(t — ti) are the convolution 
of the density of states function with the Fermi-Dirac 



distribution^!. The conductance G is equal to the cur- 
rent divided by the bias voltage. A comprehensive de- 
scription of our nume rical implementation has been pub- 
lished previousl} 122 1 49 1 . 

An exquisite many-body state called the Kondo effect 
arises when the dot level is positioned below the Fermi 
energy at sufficiently low temperatures. The net spin 
localized within the dot and the Fermi sea of electrons 
in the contacts hybridize to form a spin singlet. This 
results in a sharp resonance fixed to the Fermi levels of 
the contacts in the dot density of states. The linewidth of 
the Kondo resonance can be approximated by an energy 
scale Tk (Kondo temperature) given by 

„ (DT\* ( Tt\e Aot \ \ n7 , 

K \ ~T ) p V — ~ ) ' 

where D is a high energy cutoff equal to half bandwidth of 
the conduction electrons and Y corresponds to the value 
of coupling between the dot and the contacts Y(e) at 
e = e F . 

Our aim in this paper is to theoretically investigate 
a case in which the dot level is displaced from its equi- 
librium level abruptly by means of a gate voltage. We 
will be particularly interested in a system which has been 
studied before in the absence of any electron-phonon cou- 
pling. Its dot level is abruptly moved from edot — — 5r 
to e dot = -2T at t = where Y = Yp(e f ). We will re- 
port the instantaneous conductance right after the dot is 
moved to its final position. In the following discussion, 
we will take the phonon frequency and the bandwidth as 
ujq = 0.06r and D = 9Y respectively with T=0.4 eV. 



III. RESULTS 

We begin our analysis with the instantaneous conduc- 
tance results immediately after the dot level has been 
switched to its final position. The results shown in Fig. [2] 
correspond to three different temperatures for a constant 
nonzero electron-phonon coupling g. Instantaneous con- 
ductance results in the absence of any electron-phonon 
coupling has been reported previously^. In Fig. [2j short 
timescale corresponds to V t < 10. The conductance os- 
cillations related to the charge transf er in this timescale 
have already been analyzed in detaiP^I, therefore we 
will not elaborate on them here. The Kondo timescale 
occurs between the end of the short timescale and at- 
tainment of a plateau by the instantaneous conductance. 
Kondo resonance starts developing in this timescale. It 
takes place roughly between 10 < Y t < 60 in Fig. [2j 

The first key feature that should be mentioned re- 
garding Fig. [2] is that the steady state conductances 
(i.e. long time limit) are smaller than the case with- 
out any electron-phonon coupling for all temperatures 
as expected^. This anticipation stems from the fact 
that the electron-phonon coupling gradually quenches the 
Kondo effectH. This suppression has been explained with 
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dephasing due to electron-phonon coupling and down- 
ward shift of the energy level as a result of phonon reor- 
ganization. The second effect is an obvious consequence 
of the renormalization of the dot level. 

The second and more subtle effect is the sinusoidal os- 
cillation of the current in the Kondo timescale. This ef- 
fect is difficult to see in the main panel of Fig. [2j therefore 
in the inset, we show the magnification of the main panel 
with shifted conductance curves such that they overlap 
at the onset of oscillations. In the inset of Fig. [2] it is 
clear that the oscillation frequency is the same for all 
temperatures. Moreover, the amplitude of the oscilla- 
tions decreases as the ambient temperature increases. It 
turns out that the oscillation frequency is equal to the 
phonon frequency ujq. 



0.8 




CO I I 

50 100 150 200 

r t 

FIG. 2: Black(solid), red(dashed) and blue(dash dotted) 
curves show instantaneous conductance results as a function 
of time for T=0.0009r, T=0.0010r and T=0.0012r respec- 
tively with g=2.25 after the dot level has been switched to 
its final position. The inset is the magnified version of the 
main panel. The curves for high temperatures in the inset 
have been shifted to enable camparison of the amplitudes. 

In order to test the electron-phonon coupling strength 
dependency of the oscillation frequency, we increased g 
and performed the previous calculation at the same tem- 
peratures. The result is shown in Fig. [3] First, the steady 
state conductance is lower for all temperatures for larger 
g. This is expected due to dephasing and downward level 
shift as pointed out above. On the other hand, conduc- 
tance oscillations once more take place with a frequency 
equal to the phonon frequency ujq. We should mention 
that we increased g in Fig. [3] by increasing A and keeping 
wo constant with respect to Fig. [2] to facilitate a direct 
comparison. It must be noted that changing the value 
of loq in our calculations again yields the oscillation fre- 
quency of the conductance as luq. We checked that the 
oscillation frequency remains pinned to wo for all g val- 
ues. 

We now want to study the amplitude of oscillations for 
various electron-phonon coupling strengths and tempera- 
tures in a systematic fashion. Fig. [4] shows the behaviour 



of amplitudes as a function of electron-phonon coupling 
strength at two different temperatures. The amplitude 
is zero at g—0 regardless of the temperature but it starts 
to increase gradually until it reaches a maximum before 
g=4. It then starts decreasing and it again approaches 
zero for large g values. Meanwhile, the amplitude of os- 
cillation is always larger at lower temperature for a given 
g. We verified that this conclusion holds for other tem- 
peratures as well. 

It is of great interest to be able to provide a micro- 
scopic description for this peculiar behaviour of the tran- 
sient current. One needs to resort to the spectral function 
of the dot to this end. This can be obtained by taking 
the Fourier transform of the retarded Green's function. 
It has been shown before that the phonon sidebands are 
separated from the main Kondo peak by an integer mul- 
tiple of phonon frequency wj^^. We propose that the 
sinusoidal oscillations seen in the transient current in the 
Kondo timescale are a result of an interference between 
the main Kondo peak and its phonon sidebands. For this 
reason, the frequency of oscillation was found to be equal 
to u . For a given g parameter, the amplitude of oscil- 
lations increases as the temperature decreases because 
both the main Kondo peak and its phonon satellites are 
more developed at lower temperatures leading to stronger 
interference. 




0.2 - 



0.0 

50 100 150 200 

r t 



FIG. 3: Black(solid), red(dashed) and blue(dash dotted) 
curves show instantaneous conductance results as a function 
of time for T=0.0009r, T=0.0010r and T=0.0012r respec- 
tively with <?=3.75 after the dot level has been switched to its 
final position. 

On the other hand, sweeping electron-phonon coupling 
strength at constant T gives rise to a more complicated 
situation as seen in Fig. |4j Even though the main Kondo 
peak is developed most fully for g=0 at a given tem- 
perature T resulting in largest steady state conductance, 
no oscillation was detected in infinitesimal bias in Kondo 
timescale^ because the phonon sidebands are completely 
absent in this case. Consequently, interference is nonex- 
istent and oscillation amplitude is zero. For small g 
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values, amplitude of the oscillations gradually increases 
with g because while the main Kondo peak pinned to the 
Fermi level gets inhibited a bit, its phonon sidebands be- 
come slightly more pronounced^. This naturally leads 
to stronger interference. Note that for g <0.5, the system 
is no longer in strong coupling regime. Still, our simu- 
lations show a consistent behaviour in this limit. The 
behaviour changes for moderate values of g. The am- 
plitudes reach a maximum around <7=4 and then start 
decreasing for larger g values since destruction sets in 
for all peaks. Unsurprisingly, the amplitude goes to zero 
in large g limit where the Kondo effect disappears com- 
pletely and none of the peaks survives. 
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FIG. 4: Black(dashed) and red(dotted) curves show am- 
plitude of the second peak in conductance oscillation as 
a function of the electron-phonon coupling strength g for 
T=0.0009r and T=0.0010r respectively. The amplitudes 
have been measured with respect to the steady state (i.e. F t 
— > oo) conductance values. 

It is important to note that the oscillations take place 
exclusively with loq in infinitesimal bias and all the in- 
teger multiples of this frequency are completely absent 
as one can easily see from the Fourier transform of the 
instantaneous conductance in the long timescale shown 
in Fig. [5j This result is somewhat unexpected since the 
main Kondo peak should be interfering with other side- 
bands as well. This issue can be explained easily from an 
intuitive point of view. As one can see in Fig. [T] schemat- 
ically and confirm by taking the Fourier transform of the 
retarded Green's function, sidebands start getting drasti- 
cally smaller away from the main Kondo peak resulting in 
less interference with it. Therefore, the oscillation ampli- 
tudes associated with these peaks is negligible compared 
to the one associated with the first sideband leading to 
the absence of other frequencies. 

Before we conclude, we would like to address the ef- 
fect of finite bias on the above reported time dependent 
conductance. First, it is well known that the bias would 
split the main elastic Kondo peak into two each of which 
are pinned to the Fermi level of the contacts. These split 



Kondo resonances would induce SKP oscillations in tran- 
sient current with a frequency equal to the bias V in the 
absence of any electron-phonon coupling^. Therefore, 
the effect of bias on the above system is two fold. First, 
steady state conductance would be lower than the in- 
finitesimal bias case since the Kondo effect is gradually 
destroyed with bias. In fact, steady state investigation of 
a system with both strong electron-electron and electron- 
phonon interaction in finite bias showed that the differen- 
tial conductance exhibits enhancements when the bias is 
an integer multiple of the phonon frequency wj^. This is 
due to the fact that the inelastic phonon sidebands over- 
lap with the split Kondo peaks. Indeed, this behaviour 
is reminiscent of the time averaged conductance of an ac 
driven quantum dotP^. 

0.25i i i i 
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FIG. 5: This figure shows the Fourier transform of the in- 
stantaneous conductance in the long timescale in infinitesimal 
bias. 

In time dependent case, there are two interference ef- 
fects taking place simultaneously, namely between the 
split main Kondo peaks and between each split main 
Kondo peak and its phonon sidebands. Obviously, this 
results in two distinct oscillation frequencies. This beat- 
ing gives rise to more complicated oscillations in the tran- 
sient current as seen in Fig. [6] When we take the Fourier 
transform of the transient current, we are able to iden- 
tify these two distinct frequencies corresponding to the 
bias V and the phonon frequency ujq. This result unam- 
biguously confirms the validity of above interpretation. 
When the bias is an integer multiple of ujq, behaviour of 
the transient current is quite similar to zero bias case. In 
this case, only one frequency survives and it is equal to 

UJQ. 

IV. CONCLUSIONS 

In conclusion, the non-crossing approximation was uti- 
lized in this paper to investigate the combined effect of 
strong electron-electron and electron-phonon interaction 
on instantaneous conductance in single molecule transis- 
tor when the dot level is suddenly shifted to a position 
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FIG. 6: This figure shows the instantaneous conductance as 
a function of time for T=0.0009r with g=2.25 and voltage 
bias of V=0.03F after the dot level has been switched to its 
final position. 

where the Kondo resonance is present. Our results clearly 
showed that the instantaneous conductance displays de- 
caying sinusoidal oscillations in the long timescale in in- 
finitesimal bias. This fact set these novel oscillations 
apart from the previously predicted split Kondo peak 
oscillation!^. 



Investigation of the amplitude of these oscillations indi- 
cated that it sensitively depends on the ambient temper- 
ature and the electron-phonon coupling strength. More- 
over, the frequency of these oscillations was found to be 
precisely equal to the phonon frequency upon taking the 
Fourier transform of the instantaneous conductance in 
the long timescale. Based on these observations and in- 
vestigation of the density of states of the quantum dot, 
we proposed that the origin of this novel phenomenon can 
be attributed to the interference between the main elastic 
Kondo peak and its phonon sidebands. We also uncov- 
ered that finite source-drain bias results in more com- 
plicated fluctuation patterns in the instantaneous con- 
ductance. Two distinct frequencies corresponding to the 
phonon frequency loq and bias voltage V play role in this 
case. 

We believe the phenomenon discussed in this paper can 
be observed with present day ultrafast transient exper- 
imental techniques^ since it takes place in the longest 
timescale which typically occurs on the order of tens of 
picoseconds. Besides, the theory presented here depicts a 
fairly realistic assessment of the single molecule junctions 
as it takes into account both Coulomb and vibrational in- 
teractions that are ubiquitous for real molecules, thus we 
hope to invigorate this field by motivating a new experi- 
ment with the predictions in this paper. 
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